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^ ■ ABSTRACT 

Recently, Gamma-Ray Bursts (GRBs) were proposed to be a complementary cosmological probe 
\ to type la supernovae (SNIa). GRBs have been advocated to be standard candles since several 

O |. empirical GRB luminosity relations were proposed as distance indicators. However, there is a so- 

called circularity problem in the direct use of GRBs. Recently, a new idea to calibrate GRBs in a 
completely cosmology independent manner has been proposed, and the circularity problem can be 
> . solved. In the present work, following the method proposed by Liang et ai, we calibrate 70 GRBs 

with the Amati relation using 307 SNIa. Then, following the method proposed by Shaheloo et ai, 
we smoothly reconstruct the cosmic expansion history up to redshift z = 6.29 with the calibrated 
GRBs. We find some new features in the reconstructed results. 
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I. INTRODUCTION 

The current accelerated expansion of our universe [l| has been one of the most active fields in modern 
cosmology since its discovery in 1998 from the observations of type la supernovae (SNIa) H- Later, the 
observations of cosmic microwave background (CMB) anisotropy 3] and large-scale structure [J] confirmed 
this discovery. Although today there are already many observational methods, SNIa have been proved 
to be one of the most powerful tools to probe this mysterious phenomenon. In 1992, the famous Phillips 
relation was found @, which claimed that for nearby SNIa there exists a clear correlation between their 
intrinsic brightness at maximum light and the duration of their light curve. Since then, some empirical 
technics M, Ql H| have been developed to use the Phillips relation to make SNIa into standard candles. 
See e.g. @, [l0| for brief historical reviews. 

However, SNIa are plagued with extinction from the interstellar medium. Hence, the currently maximum 
redshift of SNIa is only about z ~ 1.755. As argued in [lflj, the observations at z > 1.7 are fairly important 
to distinguish cosmological models. On the other hand, the redshift of the last scattering surface of CMB 
is z ~ 1090. So, the observations at intermediate redshift are important. Recently, Gamma-Ray Bursts 
(GRBs) were proposed to be a complementary probe to SNIa. So far, GRBs are the most intense explosions 
observed in our universe. Their high energy photons in the gamma-ray band are almost immune to dust 
extinction. Up to now, there are many GRBs observed at 0.1 < z < 8.1, whereas the maximum redshift 
of GRBs is expected to be 10 or even larger [11| . Strictly speaking, GRBs are not standard candles, with 
radiated energies spanning several orders of magnitude. But the use of empirical luminosity correlations 
to standardize them as distance indicators has been proposed by several authors. In some sense, this is 
reminiscent of the case of the Phillips relation for SNIa. These empirical correlations include the Amati 
relation [12], those derived from it (for examples, the Ghirlanda relatio n fill, th e Yonetoku relation [ijj], 
the Liang-Zhang relation [HI], the Firmani relation [IB]), and others ll7L Il8l. fl9l [20j ] . Therefore, the 
so-called GRB cosmologyemerges recently. We refer to e.g. [13, HH, H3, l23| ] for comprehensive reviews. 

To our knowledge, in [24| GRBs were used to constrain Q m and dark energy for the first time. Similar 
works also include (25[ for examples. However, there is a so-called circularity problem in the direct use 
of GRBs [H|. In this case, to calibrate the empirical GRB luminosity relations, one need to assume 
a particular cosmological model with some model parameters a priori, mainly due to the lack of a set 
of low redshift GRBs at z < 0.1 which are cosmology independent. When one uses these "calibrated" 
GRBs (which are actually cosmological model dependent) to constrain cosmological models, the circularity 
problem occurs. To alleviate the circularity problem, some statistical methods have been proposed, such 
as the scatter method (26[, the luminosity distance method [26j], and the Bayesian method [27j ] . These 
methods were used extensively in the literature [HI, [28|, [2!| . However, they still cannot solve the circularity 
problem completely. Another method trying to avoid the circularity problem was proposed in [3(il ] , in which 
the parameters of the empirical GRB luminosity relation and the cosmological model parameters were 
fitted to the observational data simultaneously. However, for any given cosmological model, this method 
can always obtain some parameters for the cosmological model and the empirical GRB luminosity relation. 
In this sense, any cosmological model is "viable" (except for a few obviously absurd models). So, it cannot 
be used to rule out any cosmological model. Therefore, it is also an unsatisfactory method to solve the 
circularity problem completely. 

To overcome the circularity problem completely, one should calibrate GRBs in a cosmology independent 
manner. Due to the lack of a set of low redshift GRBs at z < 0.1 which are cosmology independent, it 
was proposed to calibrate the empirical GRB luminosity relation using a sufficient number of GRBs 
within a small redshift bin centered around any redshift [3l[ (see also e.g. [ll|). However, this method 
might be unrealistic, since the current sample of observed GRBs is not large enough. Recently, a new 
idea to calibrate GRBs in a completely cosmology independent manner has been proposed in |32l . [Hj| 
independently. Similar to the case of calibrating SNIa as secondary standard candles by using Cepheid 
variables which are primary standard candles, we can also calibrate GRBs as standard candles with a large 
amount of SNIa. This idea of the distance ladder was also briefly mentioned in [33]. However, in fact, [lij 
used the ACDM model rather than SNIa data themselves to calibrate GRBs and hence the calibration is 
still cosmology dependent actually. In 32], the GRB luminosity relation was calibrated with the empirical 
formula of the luminosity distance of SNIa, while the physical meaning of this empirical formula is still 
unclear and its reliability should be tested carefully. Instead, in [33j the GRB luminosity relation was 
calibrated with the interpolated distance moduli from the Hubble diagram of SNIa themselves. We refer 
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to the original papers [H,[33| for details. Using these cosmology independently calibrated GRBs, we can 
constrain the parameters of cosmological models without circularity problem. In fact, the constraints on 
ACDM model and wo — w a parameterized dark energy model have been considered in [H, HH, HH . 

In the present work, we will calibrate GRBs to reconstruct the cosmic expansion history up to redshift 
z — 6.29. The calibration method of GRB relation adopted here is the one proposed in j33|]. However, we 
will use the datasets of SNIa and GRBs which are larger than the ones used in (33[. So, more calibrated 
high redshift GRBs can be available with smaller uncertainties. In Sec. HH we use 38 GRBs at z < 1.4 
to calibrate the Amati relation while the distance moduli of these 38 low redshift GRBs are interpolated 
from 307 SNIa data [36|, [37j]. Then, the distance moduli of 32 GRBs whose z > 1.4 can be derived from 
the calibrated Amati relation. In Sec. IIII1 following the method proposed in [H, [39j], we smooth 307 
SNIa data and 32 calibrated GRBs data whose z > 1.4 to reconstruct the cosmic expansion history up to 
z = 6.29. Conclusion and discussions are briefly given in Sec. IIVI 




0.25 0.5 0.75 1 1.25 1.5 1 UJ u L 

logE Pii /(300keV) 



FIG. 1: Left panel: The Hubble diagram of 307 SNIa (black diamonds) and 38 low redshift GRBs (red stars) 
whose distance moduli are derived by using cubic interpolation. Right panel: 38 GRBs data (red stars) in the 
log_E Pi i /(300keV) — log_Bi so /erg plane. The best-fit calibration line is also plotted. See text for details. 



II. CALIBRATING GRBS WITH AMATI RELATION 



In [33[, the authors used 192 SNIa compiled by Davis et al. [40] to interpolate the distance moduli of low 
redshift GRBs. Recently, the Supernova Cosmology Project (SCP) collaboration released their latest 307 
SNIa dataset [36|, [37|]. This so-called Union compilation is the currently largest SNIa dataset. Obviously, 
a larger SNIa dataset could bring a better interpolation. On the other hand, in [33| the 69 GRBs compiled 
by Schaefer [l(j have been used. As in 10], the derived distance moduli of high redshift GRBs in [33| 
are the weighted average of all available distance moduli from five GRB luminosity relations. For each 
GRB luminosity relation, only about a dozen GRBs at z < 1.4 can be used to calibrate the corresponding- 
relation [33j]. Instead, in the present work, we calibrate GRBs only with the Amati relation, so that we 
can use a larger GRBs dataset for single GRB luminosity relation. Here, we adopt the 70 GRBs compiled 
by Amati et al. [4l[ . Notice that in [42j there is another 76 GRBs compilation which heavily overlaps the 
one of [4l|]. However, since the fluence data given in [i^] are not in the form of bolometric fluence, they 
are not convenient for our computing. As in [33j, we choose z = 1.4 to be the divide line to separate GRBs 
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into two groups, since in the 307 SNIa dataset of [36j there is only one SNIa (2003ak) whose z > 1.4. In 
the 70 GRBs compiled by Amati et al. [HI, there are 38 GRBs at z < 1.4 and 32 GRBs at z > 1.4. The 
maximum redshift of these 70 GRBs is z = 6.29 for GRB 050904. 

Several years ago, Amati et al. found the E Pt i — E[ so correlation in GRBs as E p i = K x E™ by using 
12 GRBs with known redshifts [ll|, where E p ^ = -Ep i0 bs x (1 + z) is the cosmological rest-frame spectral 
peak energy; the isotropic-equivalent radiated energy is given by 

E iso = 4 7 rd|5 bolo (l + z)- 1 , (1) 

in which Sboio is the bolometric fluence of gamma rays in the GRB at redshift z, and is the luminosity 
distance of the GRB. Later, Amati et al. have updated it i n |43| and [4l[. Up to now, some theoretical 
interpretations have been proposed for the Amati relation [21(. It might be geometrical effects due to 
the jet viewing angle with respect to a ring-shaped emission region [44j | . or with respect to a multiple 
sub-jet model structure (45J. An alternative explanation of the Amati relation is related to the dissipative 
mechanism responsible for the prompt emission [Z(| . For convenience, similar to [l(J , we can rewrite the 
Amati relation as 

log^ = A + Mog-f^-, (2) 
erg 300 kev 

where log indicates the logarithm to base 10, whereas A and b are constants to be determined. In the 
literature, the Amati relation was calibrated with the i?i so computed by assuming a ACDM cosmology 
with particular model parameters. As mentioned above, this is cosmology dependent and the circularity 
problem follows. Here, we instead use the method proposed in 33] to calibrate the Amati relation in a 
cosmology independent manner. 

As the first step, we derive the distance moduli for the 38 low redshift (z < 1.4) GRBs of 41] by using 
cubic interpolation from the 307 SNIa compiled in [3(|. We present the interpolated distance moduli \x 
of these 38 GRBs in the left panel of Fig. [TJ The corresponding error bars are also plotted. As in [33j |. 
when the cubic interpolation is used, the error of the distance modulus fj, for the GRB at redshift z can 
be calculated by 

; A? el ti J , where A l = J| ( Zj - z) / JJ (zj - Zi ) , 

in which j runs from 1 to 4 but j ^ i; on the other hand, e Mj i are the errors of the nearby SNIa whose 
redshifts are Zi. Then, by using the well-known 

fi = 51ogd L +25, (3) 

one can convert distance modulus fi into luminosity distance (in units of Mpc). From Eq. (JTJ with 
the corresponding Sboio given in Table 1 of [III], we can derive the Ei SO for these 38 GRBs at z < 1.4. 
We present them in the right panel of Fig. [T] whereas the _E Pi i for these 38 GRBs at z < 1.4 are read 
from Table 1 of [U|. Also, we present the errors for these 38 GRBs at z < 1.4, by simply using the 
error propagation. From Fig.[Tl one can see that the intrinsic scatter is dominating over the measurement 
errors. Therefore, as in [ToL [33l] . the bisector of the two ordinary least squar es [ 47T| will be used. Following 
the procedure of the bisector of the two ordinary least squares described in (47| , we find the best fit to be 

6= 1.725 and A = 52.7322, (4) 

with la uncertainties 

a h = 0.010 and cr A = 0.0065. (5) 

The best-fit calibration line Eq. © with b and A in Eq. (HJ) is also plotted in the right panel of Fig. [TJ 

Next, we extend the calibrated Amati relation to high redshift, namely z > 1.4. Since the E p ^ for the 
32 GRBs at z > 1.4 are given in Table 1 of [ll|, we can derive the £"i so from the calibrated Amati relation 
Eq. ([2|) with b and A in Eq. (jH). Then, we derive the distance moduli fi for these 32 GRBs at z > 1.4 using 
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Eqs. ((T|) and while their 5boio can be read from Table 1 of [41| . On the other hand, the propagated 
uncertainties are given by [l(| 



where 
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"log E is 
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5 ^Sbolo 
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1/2 
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'E is 



(6) 



(7) 



in which (JE iBO , s is the systematic error and it accounts the extra scatter of the luminosity relation. As 
in [HI, by requiring the x 2 /dof of the 38 points at z < 1.4 in the log-E Pi i / (300 keV) — log£?i SO /erg plane 
about the best-fit calibration line to be unity, we find that 



= 0.184. 



(8) 



We present the derived distance moduli /i with la uncertainties for these 32 GRBs at z > 1.4 in Table U 
and Fig. [21 It is worth noting that they are obtained in a completely cosmology independent manner. 




FIG. 2: The Hubble diagram of 307 SNIa (black diamonds) and 32 high redshift GRBs (red stars) whose distance 
moduli are derived by using the calibrated Amati relation. The dashed line indicates z — 1.4. 



III. SMOOTHLY RECONSTRUCTING THE COSMIC EXPANSION HISTORY 

WITH THE CALIBRATED GRBS 

Following the well-known procedure which is frequently used in the analysis of large-scale structure [48| , 
Shafieloo et al. proposed a new method in [11] to smooth noisy data directly using a Gaussian smoothing 
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function, rather than the top hat smoothing function. The iterative method calculates the lndi(z) as 



\nd L (zy n = IndiOOJUi + N(z) [^df s ( Zi ) - lnd L ( 2l )*_i] exp 



hr 



2 I l + z 



2A 2 



(9) 



where A is a constant to be given priorly, and the normalization parameter N{z) is given by 



TV(z) -1 = J^exp 



2A 2 



(10) 



In this section, following the method proposed in [38|, l39j, we smooth the 307 SNIa data and the 32 
calibrated GRBs data at z > 1.4 to reconstruct the cosmic expansion history up to z = 6.29. 
As is well known, the luminosity distance 



c(l + z) 



dz 



(ii) 



where H(z) is the Hubble parameter, and c = 2.9979 x 10 10 cm/s is the speed of light. For convenience, 
we introduce the dimensionless luminosity distance 



D L = (l + z) 



dz 

o wr 



(12) 



where E = H/H a and H = H(z = 0) = 100 h km/s/Mpc is the Hubble constant. Thus, d L = cH Q 1 D L 
and consequently 



/i = 51ogd L + 25 = 5 [log-^-+log£> L ) +25 



5 log/ + 42.3841 



In 10 



In/ + 42.3841, 



(13) 



where / = D^/h and hence InDj, = In / + Inh. On the other hand, In dL = In / + In 2997.9. Similar 
to [Isj], we can rewrite Eq. © as 



In f(z) s n = lnf(z) s n _ 1+ N(z)Y,[^f obs (^) ~ exp 

i 

where ln/ ofcs = (p° hs - 42.3841) (In 10)/5 from Eq. JEJJ). When n = 1, 
ln/(*)f = ln/(^)g + iV(z)^[ln/ obs (^)-ln/(^)g]exp 

= \nD L {z)l +N(z)J2 [ln/°" s (^) -\nD L (z t )°] exp 
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2A 2 



(14) 
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where -Dl(z)q is the dimensionless luminosity distance of the initial guess background model. The x 2 at 
any iteration is calculated as 



2 =sr [»(z i ) n - f i° bs (z i )Y 

Xn / j 
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(16) 
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FIG. 3: Computed \ 2 f° r the reconstructed results at each iteration using 307 SNIa and 32 calibrated high redshift 
GRBs at 2 > 1.4. 



The best fit is at the minimum x n - Note that the summation is over the 307 SNIa and 32 calibrated 
GRBs at z > 1.4. From Eq. {TTJl and d L = cH^D L = /c/100, we can find the Hubble parameter as 
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(17) 



Then, the deceleration parameter q(z) and the total equation-of-state parameter wtot = Ptot/ptot can be 
given by 



q(z) = (l + z) 



H'(z) 
H(z) 



1, 



Wtot{z ) = -i + -(i + z)l AJ-, 



(18) 
(19) 



where a prime denotes a derivative with respect to z. We intensively refer to the original papers [38j, |39( 
for technical details. 

As shown in [38| , the reconstructed results are not sensitive to the chosen value of A and the initial 
guess model. So, we choose A = 0.55 and use the flat ACDM model with fl m o — 0.15 as the initial guess 
model (notice that O m o = 0.15 is the best fit of (5(| and [4l|, which also fit ACDM model to GRBs data). 

For the flat ACDM model, E(z) = [O m0 (l + z) 3 + (1 - Q m0 )] 1/2 . In Fig. [3J we present the computed x 2 
for the reconstructed results at each iteration using 307 SNIa and 32 calibrated high redshift GRBs at 
z > 1.4. It is easy to see that the y 2 goes to its minimum value very fast at just fifth iteration (similar to 
the case of the first reference in |3g|). The corresponding Xmin = 322.045. In the top-left panel of Fig.|H 
we present the three reconstructed fi(z) lines with the likelihood within la, whereas the best-fit result with 
Xmin i s indicated by a blue solid line. In fact, these three reconstructed n(z) lines cannot be significantly 
distinguished. From Eqs. (fT7]l — (fT9]) . the H(z), q(z) and w to t(z) according to these three fi{z) lines are 
also plotted in the other panels of Fig. From Fig.[U we read that Hq = 68.817 km/s/Mpc, go = —0.587 
and Wtot{z = 0) = —0.724. Interestingly, while the results are familiar at relatively low redshifts which 
cover the redshift range of SNIa, the reconstructed H(z), q(z) and Wtot(z) decrease at high redshifts which 
are in the redshift range of GRBs. The Wtot crossed the phantom divide —1. 



GRB 


2 


Sbolo (10 -5 ergcm" 2 ) 


E P ,i (keV) 




050318 


1.44 


0.42 ±0.03 


115 ±25 


44.25 ±1.58 


010222 


1.48 


14.6 ±1.5 


766 ± 30 


43.97 ±1.53 


060418 


1.489 


2.3 ±0.5 


572 ±143 


45.43 ±1.60 


030328 


1.52 


6.4 ±0.6 


328 ± 55 


43.29 ±1.56 


070125 


1.547 


13.3 ±1.3 


934 ± 148 


44.47 ±1.56 


040912 


1.563 


0.21 ±0.06 


44 ±33 


43.26 ±2.08 


990123 


1.6 


35.8 ±5.8 


1724 ± 466 


44.56 ±1.61 


990510 


1.619 


2.6 ±0.4 


423 ± 42 


44.79 ±1.54 


080319C 


1.95 


1.5 ±0.3 


906 ± 272 


46.94 ±1.63 


030226 


1.98 


1.3±0.1 


289 ± 66 


44.97 ±1.59 


000926 


2.07 


2.6 ±0.6 


310 ±20 


44.38 ±1.53 


011211 


2.14 


0.5 ±0.06 


186 ± 24 


45.24 ±1.55 


071020 


2.145 


0.87±0.4 


1013 ±160 


47.81 ±1.56 


050922C 


2.198 


0.47 ±0.16 


415 ± 111 


46.83 ±1.61 


060124 


2.296 


3.4 ±0.5 


784 ± 285 


45.90 ±1.67 


021004 


2.3 


0.27 ±0.04 


266 ±117 


46.63 ±1.73 


051109A 


2.346 


0.51 ±0.05 


539 ± 200 


47.28 ±1.68 


060908 


2.43 


0.73 ±0.07 


514 ±102 


46.82 ±1.57 


050820 


2.612 


6.4 ±0.5 


1325 ± 277 


46.30 ±1.58 


030429 


2.65 


0.14 ±0.02 


128 ± 26 


46.08 ±1.57 


050603 


2.821 


3.5 ±0.2 


1333 ± 107 


47.02 ±1.53 


050401 


2.9 


1.9 ±0.4 


467 ±110 


45.75 ±1.59 


020124 


3.2 


1.2±0.1 


448 ± 148 


46.25 ±1.65 


060526 


3.21 


0.12 ±0.06 


105 ±21 


46.03 ±1.57 


030323 


3.37 


0.12 ±0.04 


270 ±113 


47.84 ±1.72 


971214 


3.42 


0.87 ±0.11 


685 ± 133 


47.45 ±1.57 


060707 


3.425 


0.23 ±0.04 


279 ± 28 


47.21 ±1.54 


060115 


3.53 


0.25 ±0.04 


285 ± 34 


47.19 ±1.54 


060206 


4.048 


0.14 ±0.03 


394 ± 46 


48.54 ±1.54 


000131 


4.5 


4.7 ±0.8 


987 ±416 


46.54 ±1.72 


060927 


5.6 


0.27 ±0.04 


475 ± 47 


48.47 ±1.54 


050904 


6.29 


2.0 ±0.2 


3178 ± 1094 


49.96 ±1.66 



TABLE I: The numerical data of 32 calibrated GRBs at z > 1.4. The first 4 columns are read from Table 1 of 
whereas the last column is derived by using the calibrated Amati relation. See text for details. 
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FIG. 4: The top-left panel is the three reconstructed lines with the likelihood within la, whereas the best-fit 
result with xLn is indicated by a blue solid line. In fact, these three reconstructed fi(z) lines cannot be significantly 
distinguished. We also plot the fi obs of 307 SNIa (black diamonds) and 32 calibrated GRBs (red stars) in the n(z) 
panel for comparison. From Eqs. (117[l — (|19p . the H(z), q(z) and Wtot(z) according to these three fi(z) lines are 
also plotted in the other panels. See text for details. 



IV. CONCLUSION AND DISCUSSIONS 



In this work, we used the cosmology independent method proposed in [33j to calibrate the GRBs and 
derived the distance moduli [i for 32 high redshift GRBs at z > 1.4, whereas the numerical results are 
given in Table HI We have used the 307 SCP Union SNIa compilation [36| and the 70 GRBs compiled 
in [4l|. Since these GRBs are calibrated in a completely cosmology independent manner, one can use the 
calibrated 32 GRBs at z > 1.4 to constrain cosmological models without circularity problem completely. 
On the other hand, following the method proposed by Shafieloo et al. [38], in this work we smoothly 
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reconstructed the cosmic expansion history up to z = 6.29 using 307 SNIa compiled in [36j and 32 
calibrated GRBs at z > 1.4. It is worth noting that this method is also model independent. 

Here are some comments on the technical details Firstly, SNIa as standard candles are affected 

by several potential problems, e.g., dust extinction, color evolution, reliability of the Phillips relation up 
to high z, uncertainty on the progenitors, possible existence of different sub-classes, etc. Thus, calibrat- 
ing GRB with SNIa propagates these uncertainties and sources of systematics also into the calibrated 
spectrum-energy correlations. In addition, even when using a cosmology independent method, in this way 
GRBs are no more an "independent" cosmological probe. In other words, by calibrating with SNIa one 
can avoid the circularity problem intrinsic in the use of GRB spectrum-energy correlations and increase 
the accuracy in the determination of spectral parameters. However, with respect to using spectrum-energy 
correlations alone, one might introduce more systematics and introduce a "circularity" with SNIa them- 
selves. Secondly, in the present work, we used the bisector of the two ordinary least squares [13] to derive 
the best-fit parameters of the calibrated E p ^ — E- lso correlation. In fact, we closely followed the method 
used by Schaefer [l(J. In this method, one changed the value of <r% iso sys by hand to force the % 2 /c?o/ of the 
38 points at z < 1.4 to be 1. One might consider that this is not a robust method. As a good alternative, 
we refer to the maximum likelihood method used, e.g., by Amati et al. [41] . 

As mentioned above, the reconstructed H(z), q(z) and Wtot(z) are interesting in some sense. To make 
these results robust, we have tried various values of A and various initial guess models. The resulted H(z), 
q(z) and Wtot(z) are still similar to the previous ones. The features in the reconstructed H(z), q(z) and 
Wtot(z) at high redshift z > 2 remain. Some remarks are in order. Firstly, one might doubt the validity 
of using GRBs as standard candles. If GRBs are not real standard candles, the reconstructed results at 
z > 1.4 are unreliable of course. In fact, the debate is not settled in the GRB community today. Some 

authors argued that GRBs cannot be used as standard candles, see m us m m, m m izi ei ei 

for examples. On the other side, some authors advocated that GRBs can be used as standard candles 

to probe cosmology, see El El El El El EE EI El EI d SEllIllMHIMIplplplMIIl 

l3lL l32l I33L l3~il l35l [54L l55U5a. I57U58L l59l. I60l.l6ll.l62l. T&3I] for instance. Today, the situation is still unclear 
and many authors choose to wait and see. Secondly, if the GRBs can be used as standard candles, since 
the sample of available GRBs at high redshifts is still not large enough, some bias might exist in the 
reconstructed H(z), q(z) and Wtot(z) and lead to the interesting features. So, before more and better 
high redshift GRBs are available, we cannot say anything conclusively. Finally, on the other side, some 
hints were found for oscillating H(z), q(z) and w to t(z) in the literature, see e.g. [H, HE HE [67], [68| and 
references therein. So, the reconstructed results obtained here might be meaningful in some sense and 
deserve further investigations. 

Obviously, due to the large scatter and the lack of a large amount of well observed GRBs, there is a l ong 
way to use GRBs extensively and reliably to probe cosmology. The Fermi Gamma-ray Space Telescope [69| 
which was launched recently might change this situation. The GRB sample will be appreciably larger, 
and the new GRBs will be better studied in some sense. In the coming Fermi era, we hope that the GRB 
cosmology could have a bright future. 
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